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Abstract 

We present an experimental demonstration of quantum optical coherence tomography (QOCT). 
The technique makes use of an entangled twin-photon light source to carry out axial optical sec- 
tioning. QOCT is compared to conventional optical coherence tomography (OCT). The immunity 
of QOCT to dispersion, as well as a factor of two enhancement in resolution, are experimentally 
demonstrated. 
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Optical coherence tomography (OCT) has become a versatile and useful biological imag- 
ing technique , particularly in ophthalmology |4(, cardiology [5J], and dermatology 
[6j]. It is an interferometric scheme that makes use of a light source of short coherence time 
(broad spectrum) []| to carry out axial sectioning of a biological specimen. Axial resolution 
is enhanced by increasing the spectral bandwidth of the source (sub-micrometer resolution 
has recently been achieved by using a light source with a bandwidth of 325 nm [8]). However, 
as the bandwidth is increased the effects of group-velocity dispersion becomes increasingly 
deleterious [9j. Various techniques have been used in attempts to counteract the effects of 
dispersion, but these require a priori knowledge of the dispersion intrinsic to the specimen 
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A quantum version o: 
recently been proposed 



OCT that makes use of an entangled twin-photon light source has 

n 

' 111 ]. A particular merit of quantum-optical coherence tomography 



(QOCT) is that it is inherently immune to dispersion by virtue of the frequency entangle- 
ment associated with the twin-photon pairs Moreover, for sources of the same 
bandwidth, the entangled nature of the twin photons provides a factor of two enhancement 
in resolution relative to OCT. 

In this letter we report the first experimental demonstration of QOCT, and show that 
the technique is indeed insensitive to group-velocity dispersion. A parallel experiment using 
conventional OCT with a source of the same bandwidth is conducted to provide a direct 
comparison of the two techniques. Using the reflections from the two surfaces of a fused- 
silica sample buried under a 10-mm-thick ZnSe window (a highly dispersive material) as 
a sample, we obtain an improvement in resolution by a factor of approximately 5. This 
improvement arises from the concatenation of two effects: dispersion cancellation and the 
factor-of-two advantage. 

We begin with a brief discussion of the principle underlying QOCT (for a comparative 
review of the theories of QOCT and OCT the reader is referred to Ref. [U(). A schematic 
of the QOCT arrangement is illustrated in Fig. The entangled twin photons may 

be conveniently generated via spontaneous parametric down-conversion (SPDC) In 
this process a monochromatic laser beam of angular frequency u p , serving as the pump, 
is sent to a second-order nonlinear optical crystal (NLC). A fraction of the pump photons 
disintegrate into pairs of downconverted photons. Both downconveted photons have the same 
polarization and central angular frequency ujq = uj p /2, corresponding to type-I degenerate 
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SPDC. We direct our attention to a noncollinear configuration, in which the photons of 
the pairs are emitted in selected different directions (modes), denoted 1 and 2. Although 
each of the emitted photons has a broad spectrum in its own right, the sum of the angular 
frequencies must always equal u p by virtue of energy conservation. 

The twin-photon source is characterized by the frequency-entangled state 



where f2 is the angular frequency deviation about the central angular frequency a>o of the 
twin-photon wave packet, C(fi) is t ne spectral probability amplitude, and the spectral dis- 
tribution S(Q) = \((Q)\ 2 is normalized such that / dflS(fl) = 1. For simplicity we assume 
S(Q) to be a symmetric function. Each photon of the pair resides in a single spatial mode, 
indicated by the subscripts 1 and 2 in Eq. (JTJ). 

The schematic illustrated in Fig. Q has, at its heart, the two-photon interferometer 
considered by Hong, Ou, and Mandel (HOM) Q]. The conventional HOM configuration 
is modified by placing the sample to be probed in one arm and an adjustable temporal 
delay (r 9 ) in the other arm. The entangled photons are directed to the two input ports 
of a symmetric beam splitter (BS). Beams 3 and 4 at its output ports are directed to two 
single-photon-counting detectors, Dx and D 2 , respectively. The coincidence rate of photons 
arriving at the two detectors, C(r q ), is recorded within a time window determined by a 
coincidence circuit (indicated by <8>). 

An experiment is conducted by sweeping the temporal delay r q and recording the in- 
terferogram C{r q ). If a mirror were to replace the sample, this would trace out a dip in 
the coincidence rate whose minimum would occur when arms 1 and 2 of the interferometer 
had equal path lengths. This dip would result from interference of the two photon-pair 
probability amplitudes, viz. reflection or transmission of both photons at the beam splitter. 

For simplicity, we neglect losses in this exposition. A weakly reflecting sample is then 
described by a transfer function H(u), characterizing the overall reflection from all structures 
that comprise the sample, at angular frequency uj: 



The quantity r(z,u) is the complex reflection coefficient from depth z and 2<f)(z,u) is the 
round-trip phase accumulated by the wave while travelling through the sample to depth z. 
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As shown previously (ll|, the coincidence rate C(r q ) is then given by 



C(r <? )cxAo-Re{A(2r 9 )}, 



(3) 



where 




dn\H(u + n)\ 2 s{n) 



(4) 



and 




(5) 



represent the constant and varying contributions, respectively. The interferogram C(r q ) 
yields useful information about the transfer function H(u) and hence about the reflectance 



The details of the QOCT experimental arrangement are shown in Fig. El For QOCT 
scans, the dotted components (mirrors Mx and M 2 , as well as detector D 3 ) are removed. 
The entangled photons, centered about Ao = 812 nm and emitted in a non-collinear config- 
uration, travel in beams 1 and 2. The photon in beam 1 travels through a temporal delay r q 
before it enters the input port of the first beam splitter, BSi. The second photon in beam 2 
goes through a second beam splitter, BS2, which ensures normal incidence onto the sample. 
The photon returned from the sample is directed to the other input port of BSi. Beams 3 
and 4, at the output of BSi, are directed to D 1 and D 2 , respectively. 

For OCT scans, the photons in beam 1 are discarded and mirrors Mi and M 2 remain in 
place (see Fig. |2j). The photons in beam 2 serve as a short-coherence-time light source. The 
reflections from the sample and mirror Mi, after recombination at beam splitter BS 2 are 
directed to detector D3 via mirror M 2 . The result is a simple Michelson interferometer, the 
standard configuration for OCT. To conduct an experiment, the temporal delay r c is swept 
and the singles rate is recorded, forming the OCT interferogram J(r c ). This arrangement 
permits a fair comparison between QOCT and OCT since both make use of a light source 
with identical spectrum. 

The initial experiment makes use of a thin fused-silica window as the sample. The transfer 
function H(u) is then given by 



where the reflectances from the front and back surfaces are |ri| 2 = |r 2 | 2 = 0.04 at normal 
incidence; L = 90 fim is the sample thickness (which is greater than the 37-/zm coherence 
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r(z,u) 11]. 



H{uj) = r 1 + r 2 e 



(6) 



length of the source), c is the speed of light in vacuum, and n pa 1.5 is the refractive index 
of the fused silica (which is taken to be independent of uj by virtue of the low dispersiveness 
of the material). Under these conditions, Eqs. (JU), ©, and (jBJ) yield 

Ao = N 2 +|r 2 | 2 (7) 

and 

A(r q ) = N 2 s(r q ) + \r 2 \ 2 s(r q - 2r d ) + 2 Re {r x r 2 * s(r q - r d ) e ^/ c }, (8) 

where s(r q ) is the Fourier transform of the source spectrum S(Q). Substituting Eqs. (J7J) and 
(JHJ) into Eq. Q yields an interferogram that contains the three varying terms in Eq. (jHJ). 

The first two terms in Eq. (jHj) are dips arising from reflections from each of the two 
surfaces. They are separated by = 2nL/c and are expected to exhibit 50% visibility since 
l r i| 2 — l r 2| 2 - The third term, which appears midway between the two dips, arises from 
interference between the probability amplitudes associated with these reflections. This term 
changes from a hump to a dip depending on the values of u p , n, L, and the arguments of r% 
and r2- 

The experimental QOCT interferogram for this sample, normalized to A , is plotted in 
Fig. 3a. The two dips, separated by the optical path length of the sample nL = 135 /im, 
exhibit 45% visibility, in close agreement with the theoretically expected value of 50%. The 
abscissa is represented in units of the scaled temporal delay cr/2, representing the physical 
displacement of the delay line, so that r stands for r q and r c alike. 

The OCT interferogram for the same sample is expected to consist of two interference- 
fringe envelopes, each with visibility calculated to be 30%, separated by r d . The experimental 
OCT interferogram for this sample, normalized to the constant background, is shown in Fig. 
3b. The centers of the envelopes, separated by nL = 135 fim, exhibit 28% visibility. 

It is apparent that the 18.5-/zm FWHM of the dips observed in QOCT provides a factor of 
2 improvement in resolution over the 37-/xm FWHM of the envelopes observed in OCT. This 
improvement, which is in accord with theory ultimately results from the entanglement 
inherent in the nonclassical light source used in QOCT. 

To demonstrate the dispersion-cancellation capability of QOCT, we bury the sample 
under a highly dispersive medium and carry out a QOCT/OCT experiment, as described 
above. The transfer function of this composite sample is then Hdi sp (u) = H (cu) e l2l3 ^ d , 
where the transfer function of the fused-silica window H{u) is given by Eq. ©, (3{ui) is the 



wave number in the dispersive medium, and d is the thickness of the dispersive medium. 
We expand (3(uq + fl) to second order in Q: j3(uo + fi) ~ /?o + + /3"fi 2 , where j3' is the 
inverse of the group velocity v at u , and /3" represents group-velocity dispersion (GVD) 

Q- 

Substituting i7 disp (u;) into Eq. (j5J), the QOCT varying term for the buried sample turns 
out to be 

Adispfo) =A(r,-2^d), (9) 

which is simply a displaced version of the result obtained for the sample in air, as provided 
in Eq. (|5j). Neither /3o nor the GVD parameter (3" appear in Eq. (^J); nor, in fact, do any 
higher even-order terms. The cancellation of GVD is an important signature of QOCT. In 
OCT 3" does not cancel and the result is a degradation of depth resolution and a reduction 

n 

of the signal-to-noise ratio [10J. 

In the experimental realization, the fused-silica window is buried beneath two cascaded 
5-mm-thick windows of highly dispersive ZnSe, slightly canted with respect to the incident 
beam to divert back-reflections, as shown at the top of Fig. 4. The GVD coefficient for 
ZnSe is /3" = 5 X 10~ 25 s 2 m~" 1 at Ao = 812 nm, which is 25 times greater than that for fused 
silica. 

Figure 4a shows the normalized QOCT interferogram for the buried sample. As predicted 
in Eq. (0), the widths of the dips remain 18.5 /mi, just as they were in the absence of the 
dispersive medium (see Fig. 3a). Thus the resolution of the QOCT scan is unaffected by 
the presence of the dispersive medium. The hump between the two dips in the QOCT 
interferogram in Fig. 4a is also unaffected by the presence of the dispersive medium. This 
would not be the case, however, if the dispersive material were between the reflecting surfaces 
rather than outside of them |ljj. On the other hand, the interference- fringe envelopes in 
the normalized OCT interferogram displayed in Fig. 4b are broadened from 37 to 92 /xm as 
a result of dispersion. 

It is also of interest to compare the visibility of the interferograms. Diverted reflection 
losses are not expected to reduce the visibility of the QOCT features, whereas this benefit 
does not accrue to OCT Q- The robustness of QOCT in this connection is evident in 
Fig. 4a, where the reduction of visibility (in comparison with Fig. 3a) arises only from 
misalignment. In Fig. 4b, on the other hand, the OCT interferogram suffers a substantially 
greater loss of visibility (in comparison with Fig. 3b), as a result of both dispersion and 



diverted losses from the four uncoated ZnSe surfaces. 

In conclusion, we have carried out proof-of-principle experiments demonstrating the suc- 
cessful operation of a new axial optical sectioning technique, quantum optical coherence 
tomography (QOCT). We have experimentally demonstrated the two principal advantages 
that stem from the frequency entanglement of the twin-photon source: dispersion cancella- 
tion and resolution doubling. 
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FIG. 1: Schematic of quantum-optical coherence tomography (QOCT). A monochromatic laser of 
angular frequency u> p pumps a nonlinear crystal (NLC), generating pairs of entangled photons. BS 
stands for beam splitter and r q is an adjustable temporal delay. Di and D2 are single-photon- 
counting detectors that feed a coincidence circuit indicated by the symbol ®. The outcome of an 
experiment is the coincidence rate C(r g ). 
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FIG. 2: Experimental arrangement for quantum/classical optical coherence tomography 
QOCT/OCT. A monochromatic Kr + -ion laser operated at X p = 406 nm pumps an 8-mm-thick 
type-I L1IO3 nonlinear crystal (NLC) after passage through a prism, P, and an aperture (not 
shown), which remove the spontaneous glow of the laser tube. BD stands for beam dump (to 
block the pump), BS for beam splitter, M for mirror, A for 2.2-mm aperture, F for long-pass filter 
with cutoff at 725 nm, and D for single-photon-counting detector (EG&G, SPCM-AQR-15). The 
quantities T q and r c represent temporal delays. For QOCT scans, the dotted components Mi, M2 
and D3 are removed, T q is swept, and the coincidence rate C(r q ) is measured within a 3.5-nsec time 
window. For OCT scans, beam 1 is discarded (beam 2 serves as the short-coherence-time light 
source), r c is swept, and the singles rate I(t c ) is recorded. 
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FIG. 3: QOCT and OCT normalized interferograms for a 90-/xm fused-silica window in air (as 
shown at top of figure). The abscissa is the scaled temporal delay cr/2, which represents dis- 
placement of the delay line (r therefore represents both r q and r c ). (a) Coincidence rate C(r q ) 
normalized to Ao (the QOCT normalized interferogram) . (b) Singles rate I(t c ) normalized to 
constant background (the normalized OCT interferogram). 
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FIG. 4: QOCT and OCT normalized interferograms for a 90-^m fused-silica window buried beneath 
two cascaded 5-mm-thick windows of highly dispersive ZnSe. As shown at the top of the figure, 
the ZnSe is slightly canted with respect to the incident beam (arrow) to divert back-reflections. 
The abscissa is the scaled temporal delay cr/2, which represents displacement of the delay line (r 
therefore represents both r q and r c ). (a) Coincidence rate C(r q ) normalized to Ao (the QOCT nor- 
malized interferogram) . (b) Singles rate I(t c ) normalized to constant background (the normalized 
OCT interferogram). 
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